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Abstract. For the first time, the combination of semi-analytic mod- 
elling of galaxy formation and A^-body simulations of cosmic structure 
formation enables us to model, at the same time, both the photomet- 
ric and the clustering properties of galaxies. Two Cold Dark Matter 
Universes provide a reasonable fit to the observed properties of galaxies, 
groups and clusters, including luminosities, colours, density and velocity 
biases. We show how the properties of galaxies and groups on small scales 
are inextricably connected with the global properties of the Universe. 



1. Introduction 

Groups of galaxies probe the intermediate scale between galaxies and clusters. 
A satisfactory theory for the formation and the evolution of these systems has 
always been difficult, because the size of groups is small enough to require the 
modelling of the properties and the internal structure of the galaxy members 
and, at the same time, large enough to require the modelling of the large scale 
structure sorrounding the groups. 

For the last twenty years, A^-body simulations have modelled groups with 
vacuum boundary conditions in order to resolve the internal kinematics of the 
galaxies (e.g. Carnevali et al. 1981; Barnes 1985; Diaferio et al. 1993; Athanas- 
soula et al. 1997). Large-scale structure simulations have provided cosmological 
boundary conditions without resolving the galaxy internal properties (e.g. Fred- 
eric 1995a, b; Nolthenius et al. 1997). Both sets of simulations have not mod- 
elled the physics of galaxy formation and evolution which provides observable 
quantities, such as luminosity and colour. However, modelling groups within a 
cosmological context with sufficient resolution is extremely important, because 
groups both trace the large-scale structure of the Universe (e.g. Ramella et al. 
1997, 1999) and are the sites of galaxy interactions (Hickson 1997; but see also 
e.g. Bettoni; Bosma; Rampazzo; Temporin, these proceedings). Galaxy inter- 
actions are connected with the formation of AGN's (Byrd, these proceedings), 
quasars (Cavaliere, these proceedings) and multiple central black holes (Valto- 
nen, these proceedings). The importance of this large-small scale connection 
has been particularly important for successfully modelling the very existence 
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of compact groups (Mamon 1986, 1989; Diaferio et al. 1994; Hernquist et al. 
1995). The enormous progress of the last few years in probing the high-redshift 
Universe, which shows a high degree of galaxy clustering (see e.g. Mazure & 
Le Fevre 1999), impels us to build a cosmologically motivated theory of galaxy 
formation. 

State-of-the-art A^-body/hydrodynamic simulations have started only re- 
cently to barely reach the resolution required to follow the formation of galaxies 
within a cosmological context, even though only some of the relevant physical 
processes are included (Pearce et al. 1999). For the time being, a more fruitful 
approach has been the combination of semi-analytic modelling of galaxy for- 
mation with A^-body simulations. The semi-analytic approach enables us to 
have full control of the relevant galaxy formation processes, which occur on par- 
sec or smaller scales, namely star formation, stellar evolution, and effects from 
supernova explosions. The A'^-body simulations yield the merging history of 
galaxies and their phase space position, which are the result of group and clus- 
ter formation occurring on megaparsec and larger scales. By combining the two 
techniques, we are able to connect the stellar population properties of galaxies 
with their clustering and kinematic properties. This approach has already been 
particularly successful (Kauffmann et al. 1999a,b; Diaferio et al. 1999). 

Here, I briefly review this technique (Section 2). I then show some remark- 
able results which the self-consistency of the approach yields automatically: the 
two component luminosity functions of groups (Section 3) and the density and 
velocity biases of cluster galaxies with differing colours (Section 4). Extracting 
mock redshift surveys from the simulation box also allows a direct comparison 
with the observed large-scale distribution of galaxies (Section 5). These results 
show the effectiveness of our approach, and also the partial success of the cos- 
mological models. 

2. The Galaxy Formation Recipe 

According to standard inflationary models (e.g. Peacock 1999). large scale struc- 
ture in the present day Universe originates from primordial density perturbations 
amplified by gravitational instability. These models also predict that most of 
the mass in the Universe is non-baryonic and therefore "dark" . Assuming a Cold 
Dark Matter (CDM) power spectrum of the perturbations (e.g. Efstathiou et al. 
1992), the clustering proceeds hierarchically: large perturbations collapse first, 
forming small dark matter halos which aggregate to form larger halos. Within 
these halos the baryonic gas cools, flows into the centre of the halo, and forms 
stars which evolve and, eventually, some of them explode as supernovae. The 
supernovac enrich the interstellar gas with metals and reheat part of the gas. 

We need two basic ingredients to implement this picture: (1) the merging 
history of dark matter halos; and (2) recipes for the relevant processes involv- 
ing the baryonic matter. Traditional semi-analytic models (e.g. Somerville & 
Primack 1999 and references therein) use Monte Carlo simulations based on the 
extended Press & Schechter (1974) formalism (e.g. Bower 1991) to derive the 
merging tree of dark matter halos. Here, instead, we derive the merger trees 
from A'^-body simulations. Moreover, we identify the galaxy harbored by the 
halo with the central dark matter particle of the halo: the galaxy and the cen- 
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tral particle share the same phase space coordinates. When two halos merge, 
the galaxy of the most massive halo "jumps" onto the central particle of the 
resulting halo. The other galaxy still retains its identity with its original dark 
matter particle, and becomes a "satellite" within the new halo. This approach 
makes the trajectory of galaxies in phase space discontinous, but it assures a 
smooth evolution of the luminosity and stellar mass of the galaxies. 

We then implement the semi-analytic technique as follows. We assume that 
the gas is in hydrostatic equilibrium within the dark matter halo. We model both 
the baryonic and the dark matter components with truncated isothermal spheres. 
The gas cools and flows on to the central particle of the halo instantaneously, 
if the cooling time is shorter than the Hubble time at that epoch, or at a given 
rate, otherwise. 

The stellar mass M* increases accordingly to 

^=aMeoidXlO^, (1) 

at -ftvir 

where Mcoid is the mass of cold gas, i^yir is the radius of the halo of mass Mvir, 
and = GMvir/i?vir is the halo circular velocity. The star formation efficiency 
a is a free parameter. 

Given the initial mass function for the stellar population (Scalo 1986), we 
can compute the number tjsN of supernovae of Type II expected per unit mass. 
We neglect supernovae of Type la. A fraction e, kept as a free parameter, of the 
kinetic energy Esn ejected by each supernova reheats part of the cold gas to 
the halo virial temperature. The mass of the reheated gas over the time At is 

reheat V J 

Satellites merge with the central galaxy on a timescale set by dynamical 
friction. We compute this timescale analytically. Our prescription does not 
allow for merging between satellites. This shortcoming is partly responsible 
for the overproduction of central bright galaxies, along with the fact that we 
do not implement a model of chemical evolution (e.g. Kauffmann &: Chariot 
1998), which affects the gas cooling rate (see Section 3). Luminosity evolution in 
different bands is computed with stellar population synthesis models (Bruzual 
& Chariot 1999). Finally, galaxy luminosities are dimmed with an empirical 
prescription for dust extinction (Cardelli et al. 1989). 

We apply our technique to two variants of a CDM Universe: a rCDM and 
a ACDM Universe. A complete description of these models and our technique 
is in Kauffmann et al. (1999a). 



3. The Luminosity Function of Galaxy Systems 

A scientific theory is required to provide testable predictions. Sometimes, it also 
automatically provides explanations for phenomena for which the theory was 
not originally constructed. 

For example, consider the luminosity function of galaxies within different 
environments. It is well known that the usual Schechter function (Schechter 
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Figure 1. Blue-band luminosity function of galaxies within sys- 
tems of different total mass: long-dashed line, M > W^'^h~^MQ; 
short-dashed line, lO^^/i-^M© < M < lO^^/j-^M©; dotted line, 
10^2/1-1^0 < M < lO^^h-^MQ-, thin solid line, M < W^'^h-^MQ. 
The bold solid line is the total luminosity function. 

1976) does not fit all galaxy samples well. Many clusters (e.g. Biviano et al. 
1995; Molinari et al. 1998; Trentham 1998), rich groups (Koranyi et al. 1998) 
and even compact groups (Hunsberger et al. 1998; see, however, e.g. de Carvalho 
et al. 1994 and Zepf et al. 1997) show a two component luminosity function: a 
roughly Gaussian distribution at the bright end, superimposed on a Schechter 
function which dominates at the faint end. Because ellipticals and SO's are 
generally more luminous than spirals and irregulars, the Gaussian component is 
mostly populated by early type galaxies, and the Schechter component by late 
type galaxies. Moreover, because this bump at the bright end peaks at different 
magnitudes depending on the system considered, when we consider a large galaxy 
sample, which includes different environments, this feature disappears, and a 
Schechter function provides a reasonably good fit (e.g. Marzke et al. 1994). 

Simulations of the formation of groups and clusters, where galaxies are re- 
solved and are allowed to merge, predict a typical two component mass function 
at later times. This result has been well known since very early A^-body sim- 
ulations (Aarseth & Fall 1980; Roos k Aarseth 1982: Cavahere et al. 1991). 
Thus, on the assumption that light is proportional to mass, galaxy merging pro- 
vides a natural explanation for the bump observed in the luminosity function of 
galaxy systems. The fact that this bump is not always present clearly contains 
information on the formation history of the galaxy system. 

Our models, designed to reproduce the global properties of galaxies and 
their large scale distribution rather than the luminosity function of individual 
clusters, also show this characteristic bump for massive systems (long and short 
dashed lines in Fig. |l]). 

Apparently, the global luminosity function (bold solid line) shows no bump. 
Note that the global luminosity function is not a Schechter function. To obtain 
the Schechter form, it is sufficient to use a dynamical friction prescription derived 
from high resolution A^-body simulations (Springel et al. 1999), which automat- 
ically includes merging between satellites (Somerville & Primack 1999), rather 
than our analytic prescription (Section 2). The agreement with the Schechter 
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Figure 2. Surface number density profiles of galaxies within clusters 
of mass M > lO^^/i~^M0 compared with observed cluster samples. 
Solid (dashed) lines show the profiles for red (blue) galaxies in the mod- 
els. Crosses and triangles are for red and blue galaxies in the CNOC 
sample (left panels) and for non-emission-line galaxies and emission- 
line galaxies for the ENACS sample (right panels). 

function can be improved further by modelling the chemical enrichment of the 
interstellar medium, which reduces the cooling efficiency at early times. 

4. Morphology-Density Relation and Velocity Bias 

The morphology-density relation has been known since the very beginning of the 
systematic investigation of galaxies (e.g. Hubble & Humason 1931). It has also 
been known that late type galaxies within clusters have larger velocity dispersion 
than early type galaxies (e.g. Moss & Dickens 1977). Modelling galaxy and 
large scale structure formation at the same time provides a self-consistent way 
of explaining these properties of cluster galaxies. 

Fig. 1^ compares our models with the surface number density profiles of 
galaxy clusters. Profiles are for blue and red galaxies (CNOC sample, Carlberg 
et al. 1997) and galaxies with emission and non-emission line spectra (ENACS, 
Biviano et al. 1997). The agreement is satisfactory. In fact, the differences 
between models and observations are simply due to the galaxy subsample def- 
inition: real galaxies have been selected according to their spectral properties 
which we do not model here. 

The models can also account for the larger velocity dispersion of blue galax- 
ies compared to red galaxies. Blue galaxies are falling on to the cluster for the 
first time and will become red later, when they have exhausted their cold gas 
reservoirs. Although this plausible explanation is well-known (e.g. Mohr et al. 
1996), our models yield this result self-consistently. Moreover, our models show 
how this result depends on the cosmological model (Diaferio et al. 1999, Diaferio 
1999). 

Our technique is therefore particularly promising. For example, when sim- 
ulations with better spatial and force resolution become available, it will be 
possible to investigate those kinematic properties of compact group members 
which, at present, seem difficult to interpret. For example, high velocity galax- 
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Figure 3. Galaxy distribution in the CfA2N catalogue projected on 
to three declination intervals. 
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Figure 4. A mock catalogue extracted from the ACDM simulation box. 



ies (Sulentic, these proceedings) may simply be due to favourable, albeit rare, 
conditions of infall from the surroundings (Diaferio et al. 1994). 



5. Mock Redshift Surveys 

Our simulation boxes have volumes that are comparable to the volume probed 
by the northern region of the Center for Astrophysics magnitude limited Red- 
shift Survey (CfA2N hereafter; Geller & Huchra 1989; Huchra et al. 1990; de 
Lapparent et al. 1991; Huchra et al. 1995; Falco et al. 1999). Therefore, we are 
able to compare our models with a large portion of the Universe by analysing 
our simulation box the same way observers analyse their data. 

We have compiled mock redshift surveys in order to have a cluster as massive 
as Coma at the same location as in the CfA2N (see Diaferio et al. 1999 for 
details). Fig. ^ shows the CfA2N and Fig. ^ shows a mock catalogue from one of 
our simulation boxes. The large scale structure in our simulation is not as sharply 
defined as in the real Universe: for example no Great Wall (Geller &: Huchra 
1989) is present in the model. Schmalzing & Diaferio (1999) have quantified 
these topological differences on large scale using Minkowski functionals. 
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Figure 5. Group abundance by velocity dispersion a. Bold lines are 
the mean number densities averaged over an ensamble of ten mock 
catalogues. Shaded areas show the 3-0" deviations. Dashed lines are 
for the CfA2N groups. Error bars on the CfA2N curves are Poisson 
3-cj deviations. Top and bottom panels are for two different model 
luminosity functions (see Diaferio et al. 1999). 

Further constraints on the model of galaxy and structure formation come 
from clustering on smaller scales, in particular from the properties of galaxy 
groups (Ramella et al. 1997, 1999; Tucker, these proceedings). 

Catalogues of groups are usually compiled using a friends-of-friends algo- 
rithm on redshift data (Huchra & Geller 1982). However, in these redshift 
surveys, we do not know the true distance to a galaxy; we only know its red- 
shift, which also includes the galaxy peculiar velocity. Therefore, the average 
properties of groups extracted from a redshift survey do not necessarily agree 
with those of groups extracted from an ideal survey in configuration space. 

A^-body simulations allow us to quantify this difference. Although previous 
analyses suggest that the differences are small (e.g. Nolthenius & White 1987; 
Moore et al. 1993; Frederic 1995a, 1995b; Nolthenius et al. 1997), these analyses 
do not include the physics of galaxy formation. 

On the other hand, we are able to compile redshift surveys that take into 
account the luminosities of galaxies, which result from their merging and stellar 
evolution histories. We find that typical parameters of the friends-of-friends 
algorithm yield kinematic properties which are in agreement with both those of 
groups selected from configuration space and those of the CfA2N groups (Fig. 
^). However, 40% of the triplets selected in redshift space are, in fact, unrelated 
galaxies which are not bound systems in configuration space. 

Some of the differences between groups in the models and in the CfA2N is 
due to the differing distribution of galaxies on large scale. An A^-body simulation 
constrained to have the same large scale structure as the local Universe (Mathis 
et al. 1999) will be valuable for further constraining the galaxy formation recipes 
using the properties of groups. 
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6. Conclusion 



Combining semi-analytic modelling with A^-body simulations is an extremely 
powerful tool for investigating galaxy and group formation within a self-consistent 
cosmological framework. A'^-body simulations with better spatial and force res- 
olution than the simulations presented here will be ideal for understanding the 
properties of compact groups, their different morphology and environment (e.g. 
de Carvalho, these proceedings), their abundance in the local Universe and at 
high redshift. 

Groups selected objectively from two-dimensional information (Prandoni et 
al. 1994), with follow-up redshift measurements (lovino, these proceedings), or 
directly selected in redshift space (Barton et al. 1996; Allam, these proceedings), 
can provide a wealth of information on both cosmology and galaxy formation 
processes. 

The physics of the intergalactic medium should also be included in our 
approach (Cavaliere et al. 1998; Cavaliere, these proceedings): both HI (e.g. 
Combes; Verdes-Montenegro; Sancisi, these proceedings) and X-ray observations 
(e.g. Mulchaey; Ponman, these proceedings) can show clear signatures of galaxy 
interactions and tell us about the dynamical state of the group (Diaferio et al. 
1995). 

To simulate the evolution of a group of galaxies, we can extract the ini- 
tial conditions and the external tidal field from a large-scale A?^-body simula- 
tion. This strategy has already been extremely successful in simulating galaxy 
clusters (Tormen et al. 1997; Springel et al. 1999). Thus, high-resolution A''- 
body/hydrodynamic simulations will enable us to predict the properties of the 
galaxies and the intergalactic medium of individual groups within a full cosmo- 
logical context. 
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